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We present a theoretical framework to accurately calculate the location of the multicrit- 
ical point in the phase diagram of spin glasses. The result shows excellent agreement with 
numerical estimates. The basic idea is a combination of the duality relation, the replica 
method, and the gauge symmetry. An additional element of the renormalization group, 
in particular in the context of hierarchical lattices, leads to impressive improvements of 
the predictions. 

1. INTRODUCTION 

Identification of the precise location of the multicritical point is an important theoretical 
challenge in the physics of spin glasses not only because of its mathematical interest but 
also for the practical purpose of reliable analyses of numerical data. The method of duality 
is a standard tool to derive th exact location of a critical point in pure ferromagnetic 
systems in two dimensions. However, the existence of randomness in spin glasses hampers 
a direct application of the duality. 

We have nevertheless developed a theory to achieve the goal by using the combination 
of the replica method, the duality applied to the replicated system, the gauge symmetry, 
and the renormalization group [1,2,3,4,5,6]. The result shows excellent agreement with 
numerical estimates. The analysis on hierarchical lattices plays a crucial role in the 
development of the theory, in particular in the introduction of the renormalization group, 
by which systematic improvements can be achieved. 

2. MULTICRITICAL POINT 

Let us consider the ±J Ising model defined by the Hamiltonian, 
H = -YsJijWj, (1) 

where is the Ising spin and denotes the quenched random coupling. The sign of Jy, 
i.e. Jij/J = Tij, follows the distribution 



Ffo) = P 8(l-r ij ) + (l-p)5{l + T ij ) 

{tf(l-T«)+<5(l + 7fc)}, (2) 



exp(K p Tij] 



2 cosh K p 
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Figure 1. A typical phase diagram of the ±J Ising model in two dimensions. The multicritical 
point (MCP) is described by a black dot and the Nishimori line is drawn dashed. 



where exp(— 2K P ) = (1 —p)/p. The multicritical point is believed to lie on the Nishimori 
line (NL) defined by K p = [3J , where (3 is the inverse temperature. See Fig. [TJ The 
restriction to the NL simplifies the problem due to the gauge symmetry [7118] . 

According to the initial theory that uses the replica method, duality and gauge sym- 
metry [Tf2"f3"f4"] . the value of p c for the multicritical point satisfies 

H( Pc ) = ~ (3) 

where H(p) is the binary entropy, — p\og 2 p — (1 — p)log 2 (l — p), for self-dual lattices. 
Equation ([31) is solved to give p c = 0.8900, which is in reasonable agreement with numerical 
estimates. The theory has also been extended to a pair of mutually dual lattices with p cl 
and p C 2 for respective multicritical points. The result is 

H(p cl ) + H(p c2 ) = l. (4) 

Hinczewski and Berker, however, found H(pi) + H(p2) = 1.0172,0.9829,0.9911 for three 
pairs of mutually dual hierarchical lattices [9]. Their values are correct to the decimal 
points shown above as one can carry out numerically exact renormalization group calcula- 
tions on hierarchical lattices. Thus Eq. (0J is a good approximation but not quite exact, 
at least for hierarchical lattices. 

3. REPLICA AND DUALITY 

Let us give a very brief summary of the theory that leads to Eqs. ([3]) and (HI). We 
generalize the usual duality argument to the n-replicated ± J Ising model. 

We define the edge Boltzmann factor (k — 0, 1, • - • , n), which represents the configuration- 
averaged Boltzmann factor for interacting spins with k antiparallel spin pairs among n 
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Figure 2. A schematic picture of the renormalization flow and the duality for the repli- 
cated ±J Ising model. 



nearest-neighbour pairs for a bond (edge). The duality gives the following relationship 
between the partition functions on the original and dual lattices with different values of 
the edge Boltzmann factors 

Z n (xo, X\, ■ ■ ■ , x n ) = Z n (xQ, x*, ■ ■ ■ , x*), (5) 

where we have assumed self duality of the lattice in that both sides share the same 
function Z n . The dual edge Boltzmann factors x* k are defined by the discrete multiple 
Fourier transforms of the original edge Boltzmann factors, which are simple combinations 
of plus and minus of the original Boltzmann factors in the case of Ising spins. 

It turns out useful to focus our attention to the principal Boltzmann factors Xq and Xq, 
which are the most important elements of the theory. Their explicit forms are 

cosh (nK + K„) *,,//- \ n 

x (K,K p ) = c l osh ^ p \ x* (K,K p ) = (V2coshK) , (6) 

where K = f3J . We extract these principal Boltzmann factors from the partition functions 
in Eq. (J5j), which amounts to measuring the energy from the all-parallel spin configuration. 
Then, using the normalized edge Boltzmann factors Uj = Xj/xo and u* = x*/x^, we have 

x (K, K p ) Nj3 z n {ui, u 2 , ■ ■ ■ , u n ) = x* (K, K p ) NB z n {u\, u* 2 ,---, «*), (7) 

where z n (ui, ■ ■ ■) and z n (u\, ■ ■ ■) are defined as Z n /xo B and Z n /(x* Q ) Na and N B is the 
number of bonds. 

We now restrict ourselves to the NL, K = K p . Figure [2] shows the relationship between 
the curves (ut(K) , u 2 (K) , ■ ■ ■ ,u n (K)) (the thin curve) and (u\(K),u* 2 (K), ■ ■ ■ ,u* n (K)) (the 
dashed curve). The arrows emanating from both curves represent the renormalization 
flows toward the fixed point C. 

The ordinary duality argument identifies the critical point under the assumption of a 
unique phase transition. We can obtain the critical point as the fixed point of the duality 
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transformation using the fact that the partition function is a single-variable function. In 
other words, the thin curve would overlap with the dashed line for such a case. 

In the present random case, on the other hand, since z n is a multivariable function, 
there is no fixed point of the duality in the strict sense which satisfies n conditions simul- 
taneously, u±(K) = ul(K),u 2 (K) = u\{K), ■ ■ ■ , u n (K) = u* n (K). This is in sharp contrast 
to the non-random Ising model. We nevertheless assume that Xq(K, K) = Xq(K, K) may 
give the precise location of the multicritical point because, when the number of variables 
of z n in Eq. (j7j) is unity (n = 1), the fixed point condition u\ = u\ implies xq = Xq. This 
relation, in the limit of n — > in the spirit of the replica method, leads to Eq. ([3]). A 
straightforward generalization to mutual dual cases gives Eq. (j4j). 

4. RENORMALIZATION GROUP ON HIERARCHICAL LATTICES 

The renormalization group provides us with an additional point of view, especially 
on hierarchical lattices. Let us remember the following features of the renormalization 
group: (i) The critical point is attracted toward the unstable fixed point, (ii) The partition 
function does not change its functional form by the renormalization on hierarchical lat- 
tices; only the values of arguments change. Therefore the renormalized system also has a 
representative point in the same space (ui(K) , u 2 (K) , ■ ■ ■ , u n (K)) as in Fig. [2j The renor- 
malization flow from the critical point p c reaches the fixed point C, (ui\ 4 , ■ ■ • , u^). 
Here the superscript means the number of renormalization steps. There is a point d c 
related to p c by the duality, which is expect to also reach the same fixed point C since p c 
and d c represent the same critical point due to Eq. (jSJ). Considering the above property 
of the renormalization flow as well as the duality, we find that the duality relates two 
trajectories of the renormalization flow from p c and from d c . The same applies to the 
whole part of both curves, thin and dashed. In other words, after a sufficient number 
of renormalization steps, the thin curve representing the original system and the dashed 
curve for the dual system both approach the common renormalized system depicted as 
the bold curve in Fig. [2} which goes through the fixed point C. 

The partition function is then expected to become a single-variable function along the 
bold curve. This fact enables us to improve the method so that the exact location of 
the multicritical point is obtained asymptotically, which can be given by Xq~* (K) = 
Xq( s ~^°°\K). If we regard xq{K) = x^(K) as the zeroth approximation for the location of 
the multicritical point, it is expected that Xfp(K) = Xq^(K) is the first approximation 
and can lead to more precise results than Xq(K) = Xq(K) does. 

Our method by the duality analysis in conjunction with the renormalization group 
indeed has given the results in excellent agreement with the exact estimations within 
numerical errors on several self-dual hierarchical lattices as summarized in Table [TJ 

5. FURTHER DEVELOPMENTS 

The above method has also been generalized to be applicable to Bravais lattices [6j. 
Let us take an example of the square lattice. Instead of the iterative renormalization, 
we consider to sum over a part of the spins, to be called a cluster, on the square lattice 
as shown in Fig. [3] to incorporate many-body effects such as frustration inherent in spin 
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p c (without RG) p c (with RG) p c (numerical) 



0. 


S900 


0. 


3920 


0.8915(6) 


0. 


3900 


0. 


3903 


0.8903(2) 


0. 


3900 


0. 


3892 


0.8892(6) 


0. 


3900 


0. 


3895 


0.8895(6) 


0. 


3900 


0. 


3891 


0.8890(6) 



Table 1 

Comparison of the methods with and without RG and numerical estimations for several 
self-dual hierarchical lattices [5J. 




cluster 1 cluster 2 



Figure 3. The basic clusters used on the square lattice. The spins marked black on the 
original lattice are traced out instead of the iterative renormalization. 

glasses. To this end, we define the principal Boltzmann factors Xq and its dual 
those with all spins surrounding the cluster in the up state. We assume that a single 
equation gives the accurate location of the multicritical point Xq (K) = x^ s \k), where 
the superscript s stands for the type of the cluster. Recent numerical investigations on the 
square lattice have given p c = 0.89081(7) [10] , p c = 0.89083(3) [IE] and p c = 0.89061(6) 
[T2"] . while the present method has estimated p c = 0.890725 by cluster 1 of Fig. [31 and 
p c = 0.890822 by cluster 2 [6j. If we deal with clusters of larger sizes, the new method 
is expected to show systematic improvements toward the exact answer from the point of 
view of renormalization. 

The method of the renormalization group is applicable also away from the NL. For 
example, the slope of the phase boundary at the pure ferromagnetic limit has been esti- 
mated to be 1/T C x dT /dp « 3.2091 • • • on the square lattice by perturbation [2j. This 
result is applicable also to any self-dual hierarchical lattices. The present method with 
the renormalization group taken into account shows that this is not the case. The result 
depends on the type of lattice, e.g. 3.2786 • ■ • and 3.4390 • ■ ■ [15J. 

6. CONCLUSION 

The hierarchical lattices provide a very effective platform to test new ideas as has been 
exemplified in the present study. Investigations are notoriously hard for spin glasses on 
finite-dimensional systems both analytically and numerically. On hierarchical lattices, 
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on the other hand, numerically exact calculations can be carried out, and, in addition, 
hierarchical lattices share many features with finite-dimensional systems in contrast to 
mean-field systems. Analytical methods can also be implemented with relative ease on 
hierarchical lattices, which leads to the significant improvements in the prediction of the 
location of the multicritical point. Hierarchical lattices will continue to play key roles in 
the studies of spin glass and other complex systems. 
We thank financial supports by the CREST, JST. 
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